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Abstract 

The occurrence of high temperature superconductivity, and the competition with magnetism, in 
stoichiometric and doped LaOFcAs and isostructural iron-oxypnictides is raising many fundamental 
questions about the electronic structure and magnetic interactions in this class of materials. There 
are now sufficient experimental data that it may be possible to identify the important issues whose 
resolution will lead to the understanding of this system. In this paper we address a number of 
the important issues. One important characteristic is the Fe-As distance (or more abstractly the 
pnictogen (Pn) height z(Pn)); we present results for the effect of 2;(Pn) on the electronic structure, 
energetics, and Fe magnetic moment. We also study LaOFeAs under pressure, and investigate the 
effects of both electron and hole doping within the virtual crystal approximation. The electric 
field gradients for all atoms in the LaOFeAs compound are presented (undoped and doped) and 
compared with available data. The observed (vr, tt, tt) magnetic order is studied and compared with 
the computationally simpler (tt, tt, 0) order which is probably a very good model in most respects. 
We investigate the crucial role of the pnictogen atom in this class, and predict the structures and 
properties of the N and Sb counterparts that have not yet been reported experimentally. At a certain 
volume a gap opens at the Fermi level in LaOFeN, separating bonding from antibonding bands 
and suggesting directions for a better simple understanding of the seemingly intricate electronic 
structure of this system. Finally, we address briefly on the possible effects of post-lanthanum rare 
earths, which have been observed to enhance the superconducting critical temperature substantially. 

PACS numbers: 
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I. BACKGROUND AND MOTIVATION 



Isostructural and isovalent LaOFeP and LaOFeAs are layered conductors, the first being 
superconducting at Tc=2.5 while the second becomes antiferromagnetically ordered at 
T^v ~ 140 K-i^ and is not superconducting. The discovery of superconductivity at 26 K in 
carrier-doped LaOFeAs^, followed by rapid improvement now up to Tc=55 in this class, 
makes these superconductors second only to the cuprates in critical temperature. Several 
dozen preprints appeared within the two months after the original publication, and many 
hundred since, making this the most active field of new materials study in recent years (since 
the discovery in MgB2, at least). 

A host of models and ideas about the "new physics" that must be operating in this class 
of compounds is appearing, pointing out the need to establish a clear underpinning of the 
basic electronic (and magnetic) structure of the system. The materials are strongly layered, 
quasi-two-dimensional in their electronic structure, by consensus. The electronic structure 
of LaOFeP was described by Lebegue,- with the electronic structure and its neighboring 
magnetic instabilities of LaOFeAs being provided by Singh and Du^. Several illuminating 
papers have appeared since, outlining various aspects of the electronic and magnetic structure 
of LaOFeAs. 

The extant electronic structure work has provided a great deal of necessary informa- 
tion, but still leaves many questions unanswered, and indeed some important questions are 
unaddressed so far. In this paper we address some of these questions more specifically. Stoi- 
chiometric LaOFeAs is AFM; then ~0.05 carriers/Fe doping of either sign destroys magnetic 
order and impressive superconductivity arises, with Tc seemingly depending little on the 
carrier concentration. Another question is: with the nonmagnetic electronic structure of 
LaOFeP and LaOFeAs being so similar, why is the former superconducting while the latter 
is (antiferro)magnetic? Surely this difference must be understood and built into bare-bones 
models, or else such models risk explaining nothing, or explaining anything. Another ques- 
tion is the effect of the structure. Unusual sensitivity to the As height 2;(As) has been noted^; 
Tc is reported to increase with applied pressure^i^^ (reduction in volume) for low values of 
doping (up to X = 0.11 in LaOi_a;Fa;FeAs, which is reported as the amount of F for optimal 
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doping); there are increases in Tc due to replacement of La with other rare earth ions, and 
the variation in size of the rare earth is often a dominant factor in the observed trends in 
their compounds. Very important also is the magnetism in these materials, as magnetism is 
a central feature in the cuprate superconductors and in correlated electron superconductors. 
Another important question is: what can be expected if other pnictide atoms can be incor- 
porated into this system: Sb (or even Bi) on the large atom side, or N on the small atom 
end. In this paper we address these questions. 

II. CRYSTAL STRUCTURE 

The members of the family of the new Fe-based superconductors crystallize in the ZrCu- 
SiAs type structure^i'^ (space group P4/nmm, Z = 2). For instance, LaOFeAs is made of 
alternating LaO and FeAs layers, as presented in Fig. [TJ The Fe and O atoms lie in planes. 




FIG. 1: (Color online) The crystal structure of LaOFeAs, showing the alternating layers of LaO 
and FeAs. 

while the As and La atoms are distributed on each side of these planes following a chess- 
board pattern. The crystal structure is fully described by the a and c lattice parameters, 
together with the internal coordinates of La and As. Experimentally, a = 4.03533 A and 
c = 8.74090 A, while z(La) = 0.14154 and z(As) = 0.6512. However to describe correctly 
the antiferromagnetic structure, a \/2a x \/2a x c cell must be used, with four Fe atoms 
per cell, as shown in full lines in Fig. [2l We will refer to this antiferromagnetic order as 
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FIG. 2: (Color online) The Qm magnetic structure of the Fe-As substructure of LaOFeAs, showing 
alternating chains of Fe spin up (red circles with black dots) and Fe spin down (blue circles with 
black crosses). The As atoms above (below) the Fe plane are represented as large (small) squares. 
The ^/2a x \/2a x c cell is represented in full lines, while the a x a x c cell is in dashed lines. 

the Qm AFM order, or equivalently as (tt, tt, 0), while the Qo AFM order corresponds to an 
antiferromagnetic order of the original cell (dashed lines in Fig. [2]) with two Fe atoms. Also, 
FM will refer to a ferromagnetic arrangement of the spins, while NM means non-magnetic. 

III. CALCULATION METHOD 

To calculate the relevant quantities, we have used density functional theory (DFT)^ii^, 
as implemented in three different electronic structure codes. The full potential local orbital 
(FPLO) code^ii^ was mainly used, while we double checked some of the calculations with 
Wien2k code^^. For most of the FPLO and LAPW calculations, the Perdew and Wang 1992 
(PW92)^ exchange-correlation (XC) functional was used, but the effect of XC functional 
was checked using also LSDA(PZ)iS, the PBE (Perdew et al 1996)^0, and another GGA 
(Perdew et al. 1992)^^ XC functionals. At each constant volume, the crystal structure 
was fully relaxed, i.e., c/a, ;z(La) and 2;(Pn) were relaxed, where Pn is the pnictogen atom. 
The errors were estimated to be within 0.5% for c/a, and 1.0% for z(La) and z(Pn). The 
relaxation was performed in the Qa/ AFM structure, with 132 irreducible k points in the 
BZ. We double checked the total energy with a finer mesh with 320 irreducible k points in 
the BZ, and the difference is very small. After relaxation, all calculations were performed 
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using dense meshes, with 320, 1027, and 637 irreducible k points in the BZ of the Qm 
AFM, Qo AFM and NM structure, respectively. In the Qm AFM structure, we used 464 
irreducible k points in the BZ to double check the result, without any noticeable difference 
in the DOS nor band structure. As for the results presented in Sect. El we used the PAW 
(projector augmented waves) method^S as implemented in the code VASP (Vienna Ab-initio 
Simulation Package }^>2i. The Perdew Burke Ernzerhof— variant of the generalized gradient 
approximation (GGA) was used for the exchange- correlation potential. A cut-off of 600 
eV was used for the plane-wave expansion of the wave function to converge the relevant 
quantities. For Brillouin zone integrations, a mesh of 9 x 9 x 7 /c-points^^ was used within 
the modified tetrahedron method^. This mesh was decreased to 9 x 9 x 3 for the cell doubled 
along the c axis. 

IV. STUDY OF LAOFEAS IN THE TETRAGONAL STRUCTURE 

LaOFeAs has a tetragonal structure (as described in Sect. [Tll) at room temperature-. 
Although it undergoes a structural phase transition at lower temperature^'i^ (see Section |V] 
), the doped (and superconducting) material LaOi-ajF^^FeAs remains in this structure down 
to low temperature, so the study of LaOFeAs in the high symmetry structure is a necessary 
step towards the understanding of the electronic structure of the whole family of compounds. 

A. Influence of XC functionals and codes on the electronic structure of LaOFeAs 

First, we studied the electronic structure of LaOFeAs in the experimental (tetragonal) 
crystal structure for different magnetic states {Qm AFM, Qq AFM, FM and NM) using 
two different codes (FPL07 and Wien2K) and different exchange-correlation functionals. 
This is necessary in view of the large number of theoretical paper8^i2ii2^i2£i22 which appeared 
recently and often contain strong disagreements. This was partly studied by Mazin et al..— 
Table [T] summarizes the results: the magnetic moment on the Fe atom together with the 
total energy differences for each magnetic state studied here. Independent of the code or 
the XC functional used, the Qm AFM state is always found to be the ground state, which 
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confirms our earlier report^. The magnetic moment for botli AFM orders are considerably 
larger than the ordered moment reported from neutron diffraction and muon spin relaxation 
experiments, while the one for the FM order is much smaller. For this last case, FPL07 
gives zero which indicates no magnetism with both PZ and PW92 XC functional; Wien2K 
gives about 0.36 with GGA and PBE and 0.13 with PW92. It appears therefore that 
the magnetic moment of Fe for the same state with different XC functionals varies by up 
to 0.5 fiB, which is unexpectedly large, although GGA is known to enhance magnetism.— 
The difference between FPL07 and Wien2K in predicting the Fe magnetic moment for each 
state may explain the total energy differences among them. Virtual doping (see subsection 
B) by 0.1 e~/Fe enhances the Fe magnetic moment in the Qm AFM state but reduces it in 
the FM state for all the XC functionals used. 

In the structural optimization (performed in the Qm state), FPL07 with PW92 (LDA) 
functional gives reasonable c/a and z(La) in good agreement with experiment, but it pre- 
dicted z(As) ~ 0.139, which is 0.011 off the experimental value, about 0.1 A in length. 
However, Wien2K with PBE(GGA) XC functional gives an optimized z(As) ~ 0.149, which 
agrees well with experimental z(As). Similar results are found in the XFe2As2 family (X=Ba, 
Sr, Ca) too. It suggests that, GGA (PBE) XC functional optimizes the FeAs-based system 
much better than LDA (PW92) XC functional. And GGA should have better performance 
in dealing with the structure (including c/a, equilibrium volume and z(As)) under pressure 
of this FeAs family. This is probably due to the layered structure of the FeAs family which 
results in large density gradient between layers, thus GGA has better description of the 
potential. But in the meantime, GGA (PBE) further overestimates the magnetic moment 
of Fe, which is already overestimated by LDA (PW92). 

B. Effect of z(As) on the electronic structure of LaOFeAs 

Then we studied how the electronic structure of LaOFeAs depends on the value of z(As). 
Table HT] shows the difference between the experimental z(As)(~ 0.150), the optimized z(As) 
(~ 0.139) and a middle value of 0.145 when using FPL07 with PW92 XC functional: decreas- 
ing z(As) (reducing the Fe-As distance) rapidly reduces the differences in energy between 
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TABLE I: Calculated magnetic moment of Fe, the amounts of total energy per Fe lie below non- 
magnetic state of FM, Qo AFM and Qm AFM states from FPL07 and Wien2K with different XC 
functionals of LaOFeAs with experimental structure. Positive A EE means lower total energy than 
NM state. 



code 


XC 


mag. mom. (/xs) 


A EE (me 


V/Fe) 






Qm 


Qo 


FM 


Qm 


Qo 


FM 


FPL07 


PW92 


1.87 


1.72 


0.00 


87.2 


24.6 







PZ 


1.70 


1.31 


0.00 


62.2 


6.9 





WIEN2k 


PW92 


1.74 


1.52 


0.13 


136.9 


78.9 







GGA 


2.09 


1.87 


0.36 


149.1 


65.2 


3.7 




PBE 


2.12 


1.91 


0.37 


158.1 


70.2 


4.5 


0.1 e~ doped 


PW92 


1.86 




0.08 


125.2 




-0.5 


0.1 e~ doped 


GGA 


2.14 




0.26 


139.7 




-0.1 


0.1 e~ doped 


PBE 


2.16 




0.27 


149.6 




2.1 



TABLE II: Calculated magnetic moment of Fe, total energy relative to the nonmagnetic (ferro- 
magnetic) states of NM/FM, Qo AFM and Qm AFM of LaOFeAs with z(As)= 0.150 (experimen- 
tal),0.145, and 0.139 (optimized) from FPL07 with PW92 XC functional. 



z(As) 


mag 


mom. (/Us) 


A EE (meV/Fe) 


Fe 3d occ.^ 




Qm 


Qo 


FM 


FM-Qm 


Qo-Qm 


maj. 


min. 


0.150 


1.87 


1.72 


0.002 


87.2 


62.6 


4.32 


2.45 


0.145 


1.70 


1.41 


0.000 


60.5 


54.0 


4.24 


2.55 


0.139 


1.48 


0.01 


0.000 


34.6 


34.6 


4.15 


2.68 



the different magnetic orderings. At z(As) = 0.145, the magnetic moments of the Qm and 
Qo states are reduced significantly in comparison with z(As) = 0.150, and the difference in 
energy has changed by around 20%, indicating important changes in the electronic structure 
upon moving the As atom. For z(As) = 0.139, the Qo AFM state has lost its moment (be- 
come the NM state), while the magnetic moment of the Qm state has decreased even more. 
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with a changing rate of 6.8 /i^/A , indicating strong magnetophonon couphng.- Therefore, 
using the experimental or optimized value for the internal coordinate of As gives quite dif- 
ferent results and might explain several of the discrepancies seen in the previously published 
works. In Figures [3] and HJ we present the corresponding band structures, total densities of 




TamlI>OS(^lai«/(-V/spinpfroclll 



FIG. 3: The bandstructure and total DOS of Qm LaOFeAs at ambient pressure computed for 
z(As)=0.150, z(As)=0.145, z(As)=0.139. 

states, and partial densities of states calculated for different values of z(As). Surprisingly, 
the band structure near Ep referred to the common Fermi level barely changes when z(As) 
decreases. Somewhat away from Ep, the bands below the Fermi level are pushed up in en- 
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ergy when z(As) is decreased, while the effect of the Fe-As distance on the bands above ep 
is less obvious, since they are pushed up or down depending on the direction of the Brillouin 
zone. For instance, along T — X and T — Z they are pushed down, so that a decrease of the 
pseudogap is expected, as shown by Fig. [31 The peaks of the DOS just above Fermi level 
move toward it when z(As) is reduced, while the DOS below the Fermi level is quite robust 
with less changes. The important decrease of the magnetic moment of Fe when the Fe-As 
distance changes is understood by looking at the Fe-3d PDOS (Fig. H]) and the last column 
of table [TTl Although the number of Fe-3d electrons remains approximately constant, the 
number of spin up electron decreases, while the number of spin down electrons is increased 
when z(As) is reduced, which overall leads to a decrease of the magnetic moment. 

C. Effect of virtual crystal doping on the electronic structure of LaOFeAs 

Since superconductivity happens only in doped LaOFeAs, it is necessary to know how 
doping will affect the underlying electronic structure and the character of each magnetic 
state. Using the experimental lattice parameters, we performed virtual crystal doping cal- 
culations on LaOFeAs using Wien2K by changing the charge of O (doping with F) and La 
(doping with Ba, but simulating doping with Sr as well), and the corresponding number of 
valence electrons. The virtual crystal method is superior to a rigid band treatment because 
the change in carrier density is calculated self-consistently in the average potential of the 
alloy. 

There is only a weak dependence of the calculated Fe magnetic moment on the electron 
doping level: 0.1 e~/Fe doping enhances it from 2.12 to 2.16 iib (see Table [T]). However, 
electron doping reduces the total energy difference (compared to NM) in both Qa/ AFM and 
FM states. The main effect of virtual crystal doping is to change the Fermi level position, in 
roughly a rigid band fashion (see the caption of Fig. [5]for more details). The band structures 
of 0.1, and 0.2 e~/Fe doped LaOFeAs in the Qm AFM phase show only small differences; 
the charge goes into states that are heavily Fe character and the small change in the Fe 3(i 
site energy with respect to that of As 4p states is minor. 

Notably, the virtual crystal approximation continues to give strong magnetic states, 
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E (eV) (Ep is indicated by the vertical lines) 



FIG. 5: Plots of undoped, 0.1 and 0.2 electron-doped LaOFeAs Qm AFM total DOS (displaced 
upward consecutively by 10 units for clarity, obtained using the virtual crystal approximation. 
Referenced to that of the undoped compound, the Fermi levels of 0.1 and 0.2 electron-doped DOS 
are shifted up by 0.20 eV and 0.26 eV, respectively. 

whereas doping is observed to degrade and finally kill magnetism and promote supercon- 
ductivity. Thus the destruction of magnetism requires some large effect not considered here, 
such as strong dynamical spin fluctuations. 

D. Electric field gradients 

We have calculated the electric field gradients (EFG) of each atom in LaOFeAs, studying 
both the effects of doping and of magnetic order. The structure used for these calculations 
is a=4.0355 A, c=8.7393 A, z(La)=0.U2, z(As)=0.650, and the PBE(GGA) XC functional 
was used in the Wien2K code. (PW92 (LDA) XC functional gives similar results and thus 
the results are not presented here.) Since the EFG is a traceless symmetric 3x3 matrix, 
only two of Vxx, '^yy, ^ zz are independent. For cubic site symmetry, the EFG vanishes, 
hence the magnitude and sign of the EFG reflects the amount and character of anisotropy 
of the charge density. For the symmetries studied here, the off-diagonal components of the 
EFG tensor for all the four atoms are zero. For the Q^/ AFM state, the \ yz component 
calculated separately for each spin for La and As is not zero, although the sum vanishes; the 
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spin decomposition gives information about the anisotropy of the spin density that is not 
available from measurements of the EFG. 

As shown in Table [Inl and Table HV] the EFGs of both Fe and As in NM and FM states are 
very similar and they are doping insensitive, except for Fe where the EFG is comparatively 
small (in tetrahedral symmetry, the EFG is identically zero). Due to the breaking of the x-y 
symmetry in the Qa/ phase, Y^x is no longer equal to V^y. In this case, the EFGs are quite 
different from those in the NM and FM states, which shows once more that the electronic 
structure in the Qm AFM order differs strongly from the ones of the NM and FM orders. 
Also, while hole doping (on the La site) and electron doping (on the O site) significantly 
change the EFG of Fe, the EFG of As is less affected. Using nuclear quadrupolar resonance 
(NQR) measurement, Grafe et ai— reported a quadrupole frequency z/q=10.9 MHz and an 
asymmetry parameter ?7=0.1 of the As EFG in LaOo.gFo.iFeAs. This observation gives Vzz 
~ 3.00 X 10^^ V/m^, which agrees reasonably well with our result of 2.6 x 10^^ V/m^ as 
shown in Table HVl in the NM state. Upon 0.1 electron or 0.1 hole doping, the EFGs are 
modified in a similar way for As but differently for Fe. 

E. Effect of pressure on the electronic structure of LaOFeAs 

Applying pressure is often used as a way to probe how the resulting effect on the electronic 
structure impacts the superconducting critical temperature and other properties. A strong 
pressure effect was shown experimentally for the members of the LaOFeAs family^i^^i^, since 
for example Tc = 43 K could be reached under pressure for LaOi_i.F^.FeAs, in case of optimal 
doping^. To begin to understand such observations, it is necessary to determine how the 
electronic structure of the parent compound LaOFeAs is changed by pressure. 

In Fig. El the magnetic moment of Fe in the Qm AFM phase versus Fe-As distance is 
presented. Two different behaviours of the magnetic moment are observed. When 2; (As) 
is varied at constant volume (zero pressure), the decrease of the magnetic moment of Fe is 
parabolic. When pressure is applied and all internal positions are optimized (hence z{As) 
changes) the change is linear until the magnetic moment drops to zero. This linear behavior 
is followed also when the As height z{As) is shifted by 0.011 to compensate for the PW92 
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TABLE III: The EFG of Fe in LaOFeAs with NM, FM and Qm AFM states at different doping 
levels from Wien2K with PBE(GGA) XC functional. The unit is lO^^ V/m^. 



Fe 






^yy 




doping 


UD 


dn 


total 


UD 


dn 


total 


NM 


undoped 


1 1 


1 1 


0.22 


1 1 


1 1 


0.22 




O.lh (La) 


0.21 


0.21 


0.42 


0.21 


0.21 


0.42 




O.le (La) 


0.01 


0.01 


0.02 


0.01 


0.01 


0.02 




O.le (0) 


0.09 


0.09 


0.18 


0.09 


0.09 


0.18 


FM 


undoped 


0.51 


-0.30 


0.21 


0.51 


-0.30 


0.21 




O.lh (La) 


0.05 


0.39 


0.44 


0.05 


0.39 


0.44 




O.le (La) 


0.31 


-0.21 


0.10 


0.31 


-0.21 


0.10 




O.le (0) 


0.31 


-0.20 


0.11 


0.31 


-0.20 


0.11 


Q,M 


undoped 


0.22 


0.03 


0.25 


-1.11 


0.54 


-0.57 




O.lh (La) 


0.60 


-1.13 


-0.43 


-1.15 


1.04 


-0.11 




O.le (La) 


-0.55 


1.00 


0.45 


-1.05 


0.24 


-0.81 




O.le (0) 


-0.54 


1.01 


0.47 


-1.07 


0.32 


-0.75 




0.2e (0) 


-0.82 


1.17 


0.35 


-1.02 


0.52 


-0.50 



(LDA) error mentioned above. Fig. [7] collects a number of results: the effect of pressure 
on the c/a ratio, the Fe-As distance, the total energy, the difference in energy between NM 
and QM states, and the magnetic moment on Fe. Under pressure, the c/a ratio, the Fe-As 
distance, and the magnetic moment of the Qm AFM state drop linearly when volume is 
reduced. The PW92(LDA) predicts an equilibrium volume of 0.925 Vq; and the total energy 
differences between NM and Q^/ AFM state gradually drops to zero at 0.78 Vq. 

The effect of pressure on the band structure is shown in Fig. [HI While the bands change 
positions under pressure, in the corresponding DOS (right panel of Fig. [8]), the first peak 
above is moved towards the Fermi level when pressure is applied, but the DOS from -0.1 
eV to is left almost unchanged by pressure. Therefore pressure should induces important 
changes in the superconducting properties of electron-doped LaOFeAs, while they should be 
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TABLE IV: The EFG of As in LaOFeAs with NM, FM and AFM states at different doping 
levels from Wien2K with PBE(GGA) XC functional. The unit is lO^^ V/m^. 



As 






^yy 




doping 


UD 


dn 


total 


UD 


dn 


total 


NM 


undoped 


fiQ 


6Q 




69 


69 


1 38 

J.. .(JO 




O.lh (La) 


0.70 


0.70 


1.40 


0.70 


0.70 


1.40 




O.le (La) 


0.65 


0.65 


1.31 


0.65 


0.65 


1.31 




O.le (0) 


0.66 


0.66 


1.32 


0.66 


0.66 


1.32 


FM 


undoped 


0.55 


0.81 


1.36 


0.55 


0.81 


1.36 




O.lh (La) 


0.58 


0.68 


1.26 


0.58 


0.68 


1.26 




O.le (La) 


0.56 


0.74 


1.30 


0.56 


0.74 


1.30 




O.le (0) 


0.58 


0.75 


1.23 


0.58 


0.75 


1.23 


Qm 


undoped 


-0.40 


-0.40 


-0.80 


0.77 


0.77 


1.54 




O.lh (La) 


-0.42 


-0.42 


-0.84 


0.68 


0.68 


1.36 




O.le (La) 


-0.41 


-0.41 


-0.82 


0.89 


0.89 


1.78 




O.le (0) 


-0.40 


-0.40 


-0.80 


0.91 


0.91 


1.82 




0.2e (0) 


-0.29 


-0.29 


-0.58 


1.03 


1.03 


2.06 



6 1.5 
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FIG. 6: Plot of the magnetic moment of Fe atom in the Qj\/ AFM state of LaOFeAs as a function 
of the Fe-As distance, both at ambient pressure and under pressure. 
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0.7 0.75 0.8 0.85 0.9 0.95 I 



V/V„ 

FIG. 7: Plot of the optimized c/a ratio, the Fe-As distances (A), the total energy of the Qm AFM 
state (eV), the total energy differences between NM and Qm AFM state (EE(NM)-EE(Qa/ AFM) 
(40 meV/Fe), the magnetic moment {fJ-s) of the Qjv/ AFM states as a function of V/Vq. 




FIG. 8: The bandstructure and total DOS of Qm LaOFeAs computed for 0.975Vo, 0.925 Vq and 
0.875 Vq. z(As) has been shifted. 

less important for hole-doped LaOFeAs. 

The Fermi surface of Qm LaOFeAs computed for different values of the volume is presented 
in Fig. [9l The first sheet is an almost perfect cylinder along the T — Z line, while the second 
sheet is made of two ellipsoidal cylinders with some bending. They appear to be very 
similar to the FS computed at ambient pressure^. The pressure has almost no effect on the 
first sheet, but it enhances the distortion of the second sheet. 
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0.975 VO 



0.925 VO 



0.875 VO 



FIG. 9: The Fermi surface of Qm LaOFeAs computed for 0.975Vo, 0.925 Vq and 0.875 Vq. z(As) 
has been shifted. 

V. INVESTIGATION OF THE STRUCTURAL DISTORTION AND OF THE 
(7r,7r,7r) MAGNETIC ORDER 

Magnetic q =(7r,7r,7r) magnetic order in stoichiometric LaOFeAs (with magnetic cell 
being ^/2 x x 2 the crystallographic cell) was reported first by de la Cruz et al.- using 
inelastic neutron scattering. Ordering takes place at Tjv = 135 — 140 K and is preceded in 
temperature by a structural distortion occurring around 155 K. These transitions have since 
been confirmed by other groups.^ A similar structural distortion was found for NdOFeAs^, 
showing that the temperature of the structural phase transition in this case is reduced by 
about 20 K in comparison with LaOFeAs. These phase transitions have been revisited^^ 
using various experimental tools (heat capacity, ultrasound spectroscopy etc.). Although 
the magnetic ordering of FeAs layers along the c axis is less likely to be crucial for the 
mechanism of superconductivity since the involved scale of energy is expected to be very weak 
is comparison with the intralayer ordering, its study is necessary to understand the complete 
system. For the same reason, and even if there are strong indications that it does not happen 
in the case of F-doped LaOFeAs, it is interesting to see whether the structural distortion 
of pure LaOFeAs can be reproduced by ab-initio calculations, and what the corresponding 
electronic structure looks like. The results we present in this section were obtained using the 
VASP code with the PBE(GGA) functional^. 
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A. (vr, vr, 0) structural order 



The structural transformational^ changes the a/2 x \pl cell (with four iron atoms; full lines 
in Fig. [2]) from tetragonal (space group PA/nmm) to orthorhombic (space group Cmma) 
or equivalently for the primitive cell (with two iron atoms; dashed lines in Fig. [2]) from 
tetragonal (space group PA/nmm) to monoclinic (space group P112/n). To simplify our 
study, the cell doubling along the c axis due to magnetic ordering is neglected for this study, 
i.e. we consider only the (vrvrO) order. We have performed a relaxation (shape of the cell as 
well as atom positions) of LaOFeAs for different volumes, the results being presented in Fig. 
[TOl The calculated equilibrium lattice parameters as well as the internal atomic positions 
are reported in Table IVl together with available experimental data. The overall agreement is 
satisfactory, the length of the a and b lattice cell vectors being slightly overestimated by our 
calculations, while the value of c is slightly underestimated. The value of \S\ (the monoclinic 
distortion angle) is overestimated by our calculations, but the very small distortion and very 
small energy difference makes this difference understandable. The important point is that 
ab-initio calculations are indeed able to reproduce the structural instability of LaOFeAs. 

As for the atom positions within the cell, the agreement is good for the positions of La, 
O, and Fe but is less satisfying for the internal position z{As) of arsenic. The difficulty 
concerning the position of As has been reported by us previously® and is related to the 
strong magnetophonon coupling that occurs in this compound. In Fig. [TOl we present the 
corresponding lattice parameters (upper plot); magnetic moment (middle plot); and internal 
coordinate of As ( zas) (lower plot), versus volume for LaOFeAs. The range of pressure 
covered goes roughly from —2.5 GPa to 2.5 GPa. By fitting the E-V data (not shown here) 
to a Birch- Murnaghan equation of state (EOS), we find LaOFeAs to have a bulk modulus 
of Bo = 75 GPa and a bulk modulus derivative Bq = 4.1. Also, from the upper plot of Fig. 
[TOl we predict that LaOFeAs is more compressible along the c axis than along the a and b 
axes, a common characteristic of layered materials. 

More important is the dependence of the magnetic moment on the volume (middle plot 
of Fig. [TOl) . This dependence has two origins: the first one is the usual dependence of the 
magnetic moments on the volume change, but in LaOFeAs, the magnetic moment on Fe 
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TABLE V: Left and middle columns: the structure parameters of LaOFeAs in its low-temperature 
phase as obtained from x-ray^ and neutron^ studies, as reported by Yildirim^^. Right column: 
results from calculations obtained after a full relaxation of a \/2 x \f2 cell with a (vrvrO) magnetic 
order, a, 6, and c are the lattice parameters, \b\ is the monoclinic distortion angle of the primitive 
cell, and La(z), As(z), 0(z), and Fe(z) are the internal coordinate of the corresponding atom. 





X-ray (120 K) 


Neutron (4K) 


Calcs. 




V2xV2 


Primitive 


Primitive 




a 


5.68262 A 


4.02806 A 


4.0275 A 


5.69 A 


b 


5.71043 A 


4.02806 A 


4.0275 A 


5.76 A 


c 


8.71964 A 


8.71964 A 


8.7262 A 


8.70 A 




0.2797° 


0.279° 


0.69° 


La(z) 


0.14171 


0.1426 


0.1418 


As(z) 


0.65129 


0.6499 


0.6451 


0(z) 





-0.0057 


0.0 


Fe(z) 


0.5 


0.5006 


0.5 



is known^ to be strongly dependent on the internal coordinate of As which changes with 
applied pressure (lower plot of Fig. [TOi) . 

The structural distortion has been addressed by Yildirim,— who approached the ques- 
tion differently and obtained different results. While our value of the Fe moment is close to 
that for the undistorted structure as would be expected, the moment reported by Yildirim is 
0.48 fiB per Fe atom. We checked carefully the possible existence of such a magnetic solution, 
but our calculations appears to be robust, with the magnetic moment of Fe being around 
2 ij,B- As a result of the different magnetic moment, his computed DOS (see Fig. 5 in Ref. 



30l ) also is different. Together with an experimental study, Nomura et al.- reported ab-initio 



calculations on LaOFeAs for both the tetragonal and orthorhombic structures, and found 
almost vanishing magnetic moments, which correspond to a non-magnetic ground-state. In 
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FIG. 10: Upper plot: Lattice parameters versus volume. Middle plot: magnetic moment versus 
volume. Lower plot: t^As versus volume. 

our case, such a state is higher in energy by about 140 meV per Fe atom for the fully relaxed 
structure, and therefore can safely be ruled out as being the true ground-state of LaOFeAs. 
The differences in calculated values that we have noted reflect an unusual sensitivity to 
details (structure, method, XC functional). 



B. (vr, vr, vr) magnetic order 

We turn now to the investigation of LaOFeAs taking into account both the true (vr, tt, tt) 
magnetic order and the structural distortion. In this case, we have used the experimental 
structural data provided by de la Cruz et al.^. As in the case of the (vr, vr, 0) order, there 
are two possible magnetically ordered states. Only one gives the (vr, vr, vr) order to be the 
ground state versus the (tt, tt, 0) order, and by only few meV per Fe atom. This small energy 
difference is near the limit of precision of our calculations, but appears to confirm the sign 
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FIG. 11: The bandstructure of LaOFeAs along high-symmetry directions in the case of a (vrvrvr) 
magnetic order for the distorted (orthorhombic) \f2 x \/2 x 2 ceU. The high symmetry points 
are defined as F (0,0,0); X (0.5,0,0); S (0.5,0.5,0); Y (0,0.5,0); and Z (0,0,0.5), in terms of 
reciprocical lattice vectors. 

of the very weak magnetic interaction along the c layers. 

The corresponding band structure is shown in Fig. [TTJ Due to the doubling of the cell 
along the c axis, there are now four bands crossing the Fermi level (see Fig. 3 of Ref. 0). 
Along S-Y there are tiny splittings around -0.25 and -0.6 eV as well as along F-X and F-Y, 
indicating the magnitude of interlayer coupling. In particular, the splitting is particularly 
large for one pair of bands just above E^? at F. Along X-S-Y, the bands nearest the Fermi 
level are hardly split at all. Overall, the band structure retains the essential features noticed 
before^, namely a pseudogap separating bonding and antibonding states over much of the 
zone, together with dispersive bands crossing the Fermi level along only one of the two 
in-plane directions (F — X, with our choice of axes). 

The total and partial densities of states are very similar to the ones in the case of a 
(ttttO) magnetic order and won't be shown here; but we notice that the rough electron/hole 
symmetry in view of the study of doped (superconducting) materials is preserved. Also, our 
calculated Fermi surface (not shown here), made of four sheets, is very similar to the one 
presented previously^ for the (vr, vr, 0) order and folded back along k^: it has two sheets along 
the Y — Z direction which are almost perfectly cylindrical, while the two other sheets are 
more distorted, but still showing a strong two-dimensional character. 
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VI. ROLE OF THE PNICTOGEN ATOM 



As mentioned at the beginning of Section I of this paper, LaOFeAs and LaOFeP are 
isostructural and isovalent, but they have quite different properties: LaOFeAs is Qm 
AFM ordered below T7v=150 K and not superconducting, while LaOFeP is a Tc=2.5 K 
superconductor^ without magnetic order. Also, they have completely different response to 
doping: either electron or hole doping will destroy the Qm AFM ordering in LaOFeAs and 
make it superconducting with Tc over 26 K- (43 K under pressure^), while in LaOFeP, dop- 
ing changes the critical temperature less significantly to only 9 K-. A deeper understanding 
of the differences of the electronic structure of these two compounds can provide insight into 
the competition between magnetic ordering and superconductivity. For similar reasons, the 
related compounds LaOFeN and LaOFeSb (although not studied experimentally yet) are 
potentially of high interest, so we also provide predictions for their electronic structure. 

Table IVII displays the experimental structure parameters for LaOFeP- and LaOFeAs^ as 
well as the predicted structure for LaOFeN and LaOFeSb after optimization (see below for 
calculation details). As a result of the increasing size of the pnictogen atom, the Fe-Pn 
length changes. In particular, the Fe-Pn distance is consistent with the sum of the covalent 
radii of Fe and Pn, which reflects the covalent bonding nature between Fe and Pn atoms in 
this family. The slight increase of the La-0 distance through the series is just a size effect 
related to the expansion of the volume . 

The values of the Fe magnetic moment for LaOFePn with FM/NM, Qo AFM, Qm AFM 
states, and their total energy differences are presented in Table IVIII Apart from LaOFeP, 
all the members of the LaOFePn family studied here have a large Fe magnetic moment in 
the Qm AFM state, the corresponding total energy being significantly lower than the ones 
corresponding to FM/NM state. 

A. LaOFeP 

LaOFeP was the first member of the iron-oxypnictide family to be reported to be 
superconductingi. The corresponding electronic structure was studied by Lebegue using ab- 
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Pn 


a (A) 


c(A) 


c/a 


z(La) 


z(Pn) 


La-0 


Fe-Pn 


Sum 


N 


3.6951 


8.0802 


2.187 


0.170 


0.109 


2.302 


2.047 


2.00 


P 


3.9636 


8.5122 


2.148 


0.149 


0.134 


2.352 


2.286 


2.31 


As 


4.0355 


8.7393 


2.166 


0.142 


0.151 


2.369 


2.411 


2.44 


Sb 


4.1626 


9.3471 


2.246 


0.127 


0.171 


2.396 


2.624 


2.62 



TABLE VI: Structural parameters of LaOFePn (Pn = N, P, As, or Sb), as obtained experimentally 
for LaOFeP— and LaOFeAs^ or from our calculations for LaOFeN and LaOFeSb. Length units are 
in A, z(La) and z(Pn) are the internal coordinate of the lanthanum atom and the pnictide atom, 
and "Sum" means the sum of Fe covalent radius and the Pn covalent radius, which is quite close 
to the calculated value in all cases. 



Pn 


mag 


mom. (/^b) 


A EE (meV/Fe) 




Qm 


Qo 


FM 


FM-Qm 


Qo-Qa/ 


N 


1.63 


0.80 


0.027 


41.0 


40.0 


P 


0.56 




0.087 


1.6 




As 


1.87 


1.72 


0.002 


87.2 


62.6 


Sb 


2.47 


2.43 


0.000 


293.8 


82.4 



TABLE VIL Calculated magnetic moment of Fe, total energy relative to the nonmagnetic (fer- 
romagnetic) states of Qo AFM, and Qm AFM states of LaOFePn from FPL07 with PW92 XC 
functional. 

initio calculations^, but considering only a non-magnetic ground-state. Since then LaOFeP 
has been studied using various experimental tools: by using photoemission^Ii^Sii^, it was 
shown that the Fe 3d electrons are itinerant, and that there is no pseudogap in LaOFeP. 
Also, magnetic measurements revealed^>^ that LaOFeP is a paramagnet, while electron-loss 
spectroscopy^^ implied a significant La-P hybridization. The absence of long-range order in 
LaOFeP was confirmed by Mossbauer spectroscopy^^ and it was proposed that LaOFeP and 
doped LaOFeAs could have different mechanisms to drive the superconductivity in these 
compounds. Also, further theoretical studies were performed^i^i^ but without studying all 
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the possible magnetic states. 

In our calculations, we find that for FM order Fe has a weak magnetic moment of about 
0.09 fiB, with a total energy very close to the NM one; this result is much like what is 
found in LaOFeAs. A remarkable difference is that the Qo AFM state cannot be obtained. 
However, we found the Qm AFM state to be the lowest in energy, but only by about 1.6 
meV/Fe, which is about two orders of magnitude less than in LaOFeAs. LaOFeP, therefore, 
presents the situation where all of the three possible magnetic states are all very close in 
energy to the nonmagnetic state, in contrast with LaOFeAs for which the Qm AFM order 
was clearly the ground state. Thus LaOFeP is surely near magnetic quantum criticality. 

The band structure of Qm AFM LaOFeP is displayed in Fig. [T^ together with total DOS 
for both Qm AFM and NM states. The band structure of Qm AFM LaOFeP is quite different 
from that of LaOFeAs with the same Qm order, with the most significant differences along 
F-X, F-Y and F-Z lines. The difference is because the breaking of the x — y symmetry is 
much smaller in the Qm AFM LaOFeP compared to LaOFeAs, because the calculated Fe 
moment is only 0.56 in LaOFeP (it is 1.87 ^b in LaOFeAs with the same calculational 
method). The corresponding DOS is also different from that of LaOFeAs: there is structure 
within the pseudogap around Fermi level in LaOFeP (See Fig. [T2l) . The difference in total 
DOS at Ei? is significant: it is only 0.2 states/eV/spin per Fe for LaOFeAs, but it is 0.6 
states/eV/spin per Fe for LaOFeP. In the NM state of LaOFeP, it is even larger with 1.6 
states/eV/spin per Fe. The DOS of Qm AFM LaOFeP is fairly flat from the Fermi level (set 
to 0.0 eV) to 0.6 eV, so that electron doping of LaOFeP will increase the Fermi level, but 
will hardly change N(Ei?) (in a rigid band picture). 

An important consequence is that there will be no expected enhancement of T^; coming 
from N(Eir) upon electron doping. In order to see a significant increase of N(Eir) in Qm AFM 
LaOFeP, an electron doping level of at least 1.2 e~/Fe is required, which seems unrealistically 
large based on the current experimental information. This conclusion remains valid in the 
case of NM LaOFeP, since apart from a peak around Fermi level, the DOS is about the same 
as for the Qm AFM state. Again, the behavior is quite different from the one of Qm AFM 
LaOFeAs: 0.1 e~/Fe doping will increase its N(Ep') by a factor of 6: from 0.2 states/eV/spin 
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FIG. 12: Plot of LaOFeP band structure in Qm AFM state and total DOS in both Qm AFM and 
NM statesg at ambient conditions with experimental lattice parameters. 



per Fe to 1.2 states/eV/spin per Fe. 




(a)FSl (b)FS2 

FIG. 13: Fermi surface of Qm AFM LaOFeP, showing the very strong differences compared to 
LaOFeAs. 

The Fermi surface of Qm AFM LaOFeP is shown in Fig. [131 Compared to the Fermi 
surface of Qm AFM LaOFeAs presented earher by Yin et al.—, the piece enclosing the F-Z 
line (containing holes) increases in size and its x — y cross section becomes more circular 
rather than elliptic. There is another piece (absent in LaOFeAs) also enclosing the F-Z 
line with the same shape but larger in size and containing electrons instead of holes. The 
two symmetric electron-type pieces of Fermi surface lying along F-Y direction in LaOFeAs 
reduces a lot in size in LaOFeP but it has two additional similar pieces lying along F-X 
direction. In LaOFeP, it has one more hole-type piece of Fermi surface surround Z point, 
which is a small cylinder. It is, understandably, quite different from the Fermi surface of 
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V/Vo 


a (A) 


c (1) 


c/a 


z(La) 


z(Sb) 


1.000 


4.092 


8.500 


2.077 


0.137 


0.165 


1.050 


4.118 


8.812 


2.140 


0.133 


0.163 


1.100 


4.141 


9.131 


2.205 


0.129 


0.161 


1.125 


4.155 


9.274 


2.232 


0.128 


0.160 


1.138 


4.163 


9.347 


2.245 


0.127 


0.160 


1.150 


4.169 


9.418 


2.259 


0.126 


0.159 



TABLE VIII: Optimized structure parameters for LaOFeSb at several volumes. The accuracy for 
c/a is within 0.3%, and within 0.8% for 2;(La) and z(Sb). 

NM LaOFeP presented earlier-. 

Therefore, while they are isostructural and significantly covalent, LaOFeP and LaOFeAs 
present quite important differences in their respective electronic structures. These differences 
must form the underpinning of any explanation of why LaOFeP is superconducting with a 
Tc which is almost electron-doping independent, while pure LaOFeAs is not superconducting 
and becomes so only upon doping. 

B. LaOFeSb 

Since the experimental crystal structure of LaOFeSb is not reported yet, we conducted 
calculations to obtain the structure. The procedure we used is the following: starting from 
the experimental volume Vq of LaOFeAs (but with As replaced by Sb), we first optimized 
c/a, z{\ja) and 2;(Sb). Then we chose a higher volume and again optimized the parameters, 
finally finding the volume that has the lowest total energy. Using this scheme, the optimized 
volume is 1.046 Vq while for LaOFeAs the equilibrium volume is about 0.919 Vq. Assuming 
that PW92 overbinds equally for LaOFeSb as for LaOFeAs, the experimental equilibrium 
volume for LaOFeSb should be 1.046/0.919=1.138 Vq. Therefore, we performed calculations 
for a range of volume from V = Vq to V = 1.150 Vq, the corresponding structural parameters 
being presented in Table IVIIII 
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FIG. 14: Plot of Qm AFM LaOFeSb band structure and total DOS at 1.138 Vq with both optimized 
and shifted z(Sb). 

Since for LaOFeAs in the Qm AFM phase PW92 underestimated 2;(As) by 0.011 at its 
experimental volume, we corrected z(Sb) by adding 0.011 to the optimized 2;(Sb) (we refer 
to this position at the ""shifted z(Sb)"). Both for the NM and Qm AFM case, there are 
very small differences near Ep between the optimized z{As) and shifted z{As) in the band 
structure and DOS, as seen in Fig[T31 However, shifting z{Sh) induces important changes 
in the energy differences between NM and Qm AFM states, as shown in Table HXl Also, the 
magnetic moment of Fe, and the energy differences among NM/FM, Qq AFM and Qm AFM 
are strongly dependent on the volume. With decreasing volume, the difference in energy 
between the different magnetic states decreases quickly. 

At 1.138 Vo, the inferred equilibrium volume of LaOFeSb, the properties of NM/FM, Qo 
AFM, and Qm AFM are very similar to the ones of LaOFeAs at its experimental volume. 
Thus from these results we expect that doped LaOFeSb should have similar properties (viz, 
value of Tc) as LaOFeAs. 

C. LaOFeN 

The structure of LaOFeN is also not reported experimentally. In order to obtain it, the 
same procedure as for LaOFeSb was used. The lowest total energy is at 0.762 Vq' (here Vq' 
is the experimental volume of LaOFeP.). Again assuming PW92 makes a similar error as it 
makes in LaOFeAs, we estimate its equilibrium volume to be close to 0.825 Vq'. At 0.825 
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V/Vo 


mag 


mom. (/ifi) 


A EE (meV/Fe) 




Qm 


Qo 


FM 


FM-Qa/ 


Qq-Qm 


1.000 


1.58 


1.12 


0.36 


60.1 


60.1 


1.050 


1.87 


1.74 


0.44 


95.6 


68.0 


1.100 


2.09 


2.00 


0.00 


147.6 


70.7 


1.125 


2.17 


2.10 


0.00 


172.6 


71.8 


1.138 


2.23 


2.16 


0.00 


190.6 


72.5 


1.150 


2.26 


2.19 


0.00 


199.0 


72.0 


1.050 


2.17 


2.08 


0.72 


158.1 


78.0 


1.100 


2,35 


2.00 


0.00 


223.8 


80.8 


1.125 


2.42 


2.37 


0.00 


271.6 


81.8 


1.138 


2.47 


2.42 


0.00 


293.8 


82.4 


1.150 


2.49 


2.45 


0.00 


287.6 


82.1 



TABLE IX: Calculated magnetic moment of Fe, total energy relative to the nonmagnetic (ferro- 
magnetic) states of Qo AFM and Qa/ AFM with the optimized structure of LaOFeSb at several 
volumes from FPL07 with PW92 XC functional. Upper part: z(Sb) is optimized. Lower part: 
z(Sb) is optimized and shifted. 

Vo' and for larger volume, the total energy of the Qm AFM state is well below that of the 
FM/NM state (see Table IX]). Therefore, LaOFeN, if it exists, should be in the Qm AFM 
ordered state at low temperature, which is similar to LaOFeAs and LaOFeSb. 

Compared to the other LaOFePn compounds, LaOFeN is even closer to being a semimetal 
when the volume is equal to 0.825 Vq', and it becomes a small gap insulator at 0.850 Vq' 
and a higher carrier density metal at 0.800 Vq' (see Fig. [T5l) . The DOS for 0.825 Vq' shows 
a pseudogap around Ep, but the DOS is somewhat less flat than it is for LaOFeAs. 

When LaOFeN is calculated to be insulating (for volumes larger than 0.825 Vq'), the gap 
can be taken to define a distinction between bonding (occupied) and antibonding (unoccu- 
pied) states. The appearance of this gap in LaOFeN is quite surprising: although there is 
clear separation of valence and conduction bands over most of the zones for LaOFeAs, there 

27 



V/Vo' 


mag 


mom. {/j-b) 


A EE (meV/Fe) 




Qm 


Qo 


NM/FM 


NM/FM-Qa/ 




0.900 


2.21 


1.69 


1.64 


209.8 


135.9 


0.875 


2.06 


1.51 


0.03 


114.9 


99.2 


0.850 


1.88 


1.14 


0.03 


74.3 


68.1 


0.825 


1.63 


0.80 


0.03 


41.0 


40.0 


0.800 


1.26 




0.00 


18.4 




0.787 


1.08 




0.00 


11.3 




0.775 


0.90 




0.00 


7.0 




0.762 


0.00 




0.00 


1.3 




0.750 


0.00 




0.00 


1.4 




0.725 


0.00 




0.00 


1.2 




0.700 


0.00 




0.00 


0.9 





TABLE X: Calculated magnetic moment of Fe in LaOFeN, total energy relative to the nonmag- 
netic (ferromagnetic) states of Qo AFM and Qm AFM with the optimized structure at several 
volumes, but shifted z(N) up by 0.011, as a compensation PW92 does to LaOFeAs, where PW92 
underestimates z(As) by 0.011. 




FIG. 15: Plot of LaOFeN Qm AFM band structure and total DOS at 0.850Vo', 0.825 Vq' and 
0.800 Vo' with shifted z(N). 
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is no way to ascribe the small FSs to simple overlapping valence and conduction bands: in 
LaOFeAs and LaOFeSb, the bonding and antibonding bands are never completely separated 
from each other. In LaOFeN this separation finally becomes apparent, as an actual bandgap 
does appear. 

VII. ROLE OF THE RARE EARTH ATOM IN REOFEAS 

After LaOFeAs was discovered, after appropriate variation of the carrier concentration, to 
be superconducting at 26 K, much substitution on the rare earth (R) site has been done, with 
impressive increases in the critical temperature. Since all are evidently trivalent and donate 
three valence electrons to the FeAs layer, it becomes important to uncover the infiuence of 
the R atom: is it some aspect of the chemistry, which does differ among the rare earths? is 
it an effect of size? or can there be some other subtle effect? 

Table IXII is a collection of the lattice constants a and c, volume V of the primitive cell, 
Tc onset of ROFeAs reported from experiment.— '^'^'^ Both lattice constants, hence the 
volume, decrease monotonically as the atomic number increases, but increases only from 
La to Gd, whereupon drops for heavier rare earths. Since we have found that small details 
affect the electronic and magnetic structure - especially 2;(As) - it is reasonable to assess the 
size effect. We have performed calculations on Ce, Nd and Gd, using LSDA+U with U=7.0 
eV and J=l eV applied to the R atom to occupy the 4/ shell appropriately and keep the 
4/ states away from the Fermi level. Our results indicate that all have very similar DOS 
and band structure with LaOFeAs. To investigate further, we checked GdOFeAs using the 
crystal structure of LaOFeAs. The resulting band structure and DOS are almost identical 
to the original results for Gd, thus there seems to be no appreciable effect of the differing 
chemistries of Gd and La. This negative result supports the idea that the size difference 
may be dominant, though seemingly small. The difference in size (hence a, c, and the 
internal coordinates) influences not only the band structure and DOS, but also the magnetic 
properties. Fixed spin moment calculations in the FM state gives the lowest total energy at 
0.2 /ifi/Fe in LaOFeAs, and 0.5 /i^/Fe in both GdOFeAs and La-replaced GdOFeAs. 
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TABLE XI: Collection of the lattice constants a (A) and c(A), volume V (A^ of the primitive cell, 
Tc onset s (onset, middle, and zero, in K) of ReOFeAs reported from experiments. 





7 

Li 






V (^A j 




T a 

Lj<x 


/ 


4t.Uoo 


o. / oy 


1 zL9 1 zL 






Oo 




o.DOZ 


1 '^S 9Q 


40.0 


Pr 


59 


3.985 


8.595 


136.49 


51.3 


Nd 


60 


3.965 


8.572 


134.76 


53.5 


Sm 


62 


3.933 


8.495 


131.40 


55.0 


Gd 


64 


3.915 


8.447 


129.47 


56.3 


Tb 


65 


3.899 


8.403 


127.74 


52 


Dy 


66 


3.843 


8.284 


122.30 


45.3 



VIII. SUMMARY 

We have investigated in some detail the electronic structure and magnetic properties of 
the LaOFeAs class of novel superconductors using ab-initio methods. The effects of the 
Fe-As distance, of doping, and of pressure, as well as calculations of the EFGs have been 
reported. It was found that (approximate) electron-hole symmetry versus doping, and strong 
magnetophonon coupling are primary characteristics of the LaOFeAs system, and are two of 
the ingredients that need to be understood to proceed toward the discovery the mechanism of 
superconducting pairing. We studied effects of the structural distortion and of the (7r,7r, tt) 
magnetic order, finding that experiments can be reproduced fairly well by our calculations. 
Finally, the related materials LaOFeP, LaOFeSb, and LaOFeN were investigated and their 
properties compared to those of LaOFeAs. From these comparisons, it appears that LaOFeP 
is significantly different from the other materials studied here; this difference might explain 
why, at stoichiometry, LaOFeP is superconducting while LaOFeAs is antifcrromagnetic. Also, 
in view of their similarities with LaOFeAs, either pure or doped LaOFeSb and LaOFeN are 
potential candidates as superconductors. 
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